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We show that the Deterministic Lattice Gas (DLG) (Phys. Rev. Lett. 64, 3103 (1990)) model of 
Self-Organized Criticality (SOC) despite of its deterministic micro dynamics belongs to the Manna 
universality class of absorbing state phase transitions. This finding is consistent with our observation 
that an effective stochastic term is generated in the DLG at large length scales, whereby the macro 
dynamics of the DLG appears closer to the other stochastic SOC models of the Manna class. 



It is well known that Self-Organized Criticality (SOC) 
was introduced by Bak, Tang and Wiesenfeld (BTW) as 
an attempt to explain the widespread occurrence of 
temporal fluctuations and fractal spatial structure 
It was, however, soon realized that the model used by 
BTW did not contain a 1// speetrumjl, i]. The De- 
terministic Lattice Gas (DLG) was introduced in an ef- 
fort to find another deterministic model, which could 
be used as a proof of existence for the SOC scenario. 
The model being deterministic is difficult to analyze an- 
alytically and one has to rely on a combination of sim- 
ulations and non-rigorously justified effective analytic 
investigations Numerical simulations found that the 
density fiuetuations in the DLG did exhibit 1// fluctua- 
tions (and that dissipation did take place on a fractal). 
This was found to be consistent with analysis of (non- 
linear) diffusion equations for which one assumes the ab- 
sence of any bulk noise term[3|. 

Until now it has been unclear how the DLG related 
to the universality classes of SOC[^ Q and how the ex- 
istence of an absorbing state at low density influences 
the behavior of the model. In this letter we flrst present 
simulations which demonstrate that the density fluctua- 
tions at elevated densities for large systems have a power 
spectrum S{f) oc 1//'* with /i = 3/2, which is consistent 
with a description in terms of diffusion equation with a 
bulk noise term. Next we present compelling evidence 
that the transition to an absorbing state is of the Manna 
universality class and use scaling relations to determine 
the value of the power spectrum exponent /i for densities 
near the absorbing state phase transition. 

Model - In the DLG particles interact through a 
nearest-neighbor repulsive central unit force and double 
occupancy is not permitted. All particles are updated 
simultaneously by moving the particles deterministically 
to neighbor sites according to the vector sum of the forces 
they are subject to. If two particles want to move to the 
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same site, the particle subject to the strongest force is 
moved, while in the case of equal forces no particle is 
displaced. Periodic boundary conditions are considered 
and we consider the model in dimension d — 2. 

The number of particles N(t) in a sub-volume of the 
lattice exhibits interesting temporal fluctuations. We de- 
termine the power spectrum S{f) oi N{t) from the square 
of the absolute value of the Fourier transform. Successive 
time sequences are averaged in order to achieve sufficient 
statistics: 



Sif) = ^ 
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The angular brackets denote averaging over many differ- 
ent time series. 

The spectrum S{f) was shownQ to satisfy S{f) ^ 
1//^ at the particle density p = 0.3, while previous works 
showed that the same result is obtained in a wide range of 
densities with a drive at the boundary Q. In those papers 
the maximum linear system size considered is L = 128, 
we will see that different behaviors are observed for larger 
systems. In Fig. [T] we show how the effective exponent 
/i changes from = 1 to /x = 3/2 as L increases. It is 
interesting to relate this change in the temporal fluctua- 
tions to the appropriate macroscopic equation of motion 
for the particle density p{r,t). 

Simulations show that although particles move deter- 
ministically in the DLG model they do behave like ran- 
dom walkers in the sense that their square displacement 
is linear in time. This suggests that p(r, t) evolves ac- 
cording to a diffusion equation and integrating the diffu- 
sion equation allows us to extract the temporal evolution 
of N{t) = Jy drp{r,t), where V denotes the measuring 
sub-volume. Assume that p{r, t) satisfies the diffusion 
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It has been shown0,|3l that the power spectrum of N{t) 
depends on the details of the diffusion equation. If the 
equation is driven by white noise at the boundary of 
V and the noise term ^ is absent the power spectrum 
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FIG. 1: Scaling behavior of the spectrum <S'(/) oc /"'' of the 
total number of particles N{t) in the DLG for increasing linear 
sizes of the lattice L. A crossover from — 1 for small L to 
fi = 1.5 for large L is observed. Particle density p = 0.5. 
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FIG. 2: Scaling behavior of the spectrum <S'(/) oc f~'^ of the 
total number of particles N{t) in the DLG for different par- 
ticle densities p. A crossover from p ~ 1.8 at p ~ pc — 0.245 
to p — 1.5 at p ^ pc is observed. S{f) has been multiplied 
by different constants for different densities p to visualize the 
scaling exponents properly. Lattice size L = 1000. 



exponent derived from Eq. (I2|) is /i = 1. If in contrast 
a conservative bulk noise term is included in Eq. ^ 
one obtains instead pL = 3/2. Including non-linearities 
in Eq. ^ will not influence the behavior of p,, see [1]. 
It is this observation that allows us to conclude that a 
conservative noise term is generated in the macroscopic 
Langevin description when the sub-volume V becomes 
sufficiently large. 

S{f) scaling for large lattices - We now consider the 
DLG for much larger lattices than previously studied and 
our aim is to understand the fluctuations spectrum as 
function of density. We plot S{f) in Figure [2] and find 
that the value of /i depends on the density: 



S{f) 
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where Pc is the critical density of the absorbing phase 
transition (APT in the following). 

The observed exponent at high densities is the same 
of a gas of random walkers , while at low densities the 
power spectrum scaling exponent is determined by the 
critical properties of the DLG at the APT, as we shall 
discuss in the following. 

Absorbing Phase Transition - At very low densities 
the DLG enters configurations in which all particles are 
far away from each other and, due to the short range 
interactions, the particles become unable to move, i.e. 
the dynamics is frozen. Near the transition the density 
of active particles pa behaves according to Pa{Sp) ^ 5p^ 
for p > Pc and Pa{Sp) = for p < pc with Sp = p — pc- 
We determine the critical density pc from a log-log plot 
of Pa versus Sp and obtain pc = 0.24500(2). A regression 
analysis yields the value of the order parameter exponent 



P = 0.634(2). 

We find that close to the critical point the fluctuations 
of the order parameter Apa = L? [(p^) — (pa)^] ^ Sp^^ 
with i = 0.40(2). 

To circumvent finite size effects we make use of the 
implementation of the external field for APTs developed 
in [T^]: at each time step, after performing the DLG 
update, we choose randomly hi? particles on the lattice 
and move each of them to one of its empty neighbors. 
The external field prevents the system from falling into 
an absorbing state. We assume that the order parameter 
and its fluctuations satisfy pa{Sp,h) = R{5pX,h\'') 
and Apa{6p,h) = X^' b{5p\hX'^), where A > 0. From 
this scaling ansatz we conclude that Pa{Sp = 0,h) ^ h^/°' 
and Apa{5p = 0, ft-) h^'' . Our simulations lead to 
cr = 2.19(1) [m and 7' = 0.37(1). 

be p r 

see Fig. [3] for Ap and [l^l for pa and more details. 

Finite Size Scaling - We assume the fol- 
lowing finite size scaling form pa{Sp,h,L) = 
X-^Rpbc{SpX,hX'',L\-''^) and Apa{6p,h,L) 
X'^ Dpbc{SpX,hX'^ jLX^"^), where the exponent 
describes the divergence of the spatial correlation length, 
i.e., oc 6p~'^^. The universal scaling functions depend 
on the particular choice of boundary conditions, although 
in the thermodynamic limit Rpbdx, y, 00) = R{x, y) and 
Dpbc{x,y, 00) = D{x,y). 

Following we consider the fourth order cumulant 



Data collapses can be p roduced choosing X = h 



Q, which is defined as: Q = 1 



For non- 



vanishing order-parameter the cumulant tends to Q = 
2/3 in the thermodynamic limit. One expects that 
the cumulant obeys the scaling form: Q{Sp, h, L) = 
Qpbc{6pX,hX'',LX-''^). Choosing A = L^/"^ at 6p = 
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FIG. 3: Data collapse of Apa- In the inset data before the 
collapse are shown. Lattice linear size L = 1000. 



we obtain the following equation; 

Q(0>,i) = Qpfc,(0,/iL'^/'^^,l) 



(4) 



which enables us to determine v±^ through a data collapse 
by plotting Q{0,h,l) against hL'^/"^ as in Fig. H) Best 
results are obtained for i'±^ = 0.83(5). 
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FIG. 4: (lin-log) Data collapse of Q{0, h, L) for the determi- 
nation of v±. In the inset data before the collapse are shown. 



Dynamical scaling - Starting simulations of the DLG 
from a random distribution of particles above the crit- 
ical point pc the density of active sites pa decreases in 
time and tends to its steady state value. At the critical 
point p ~ Pc the order parameter decays algebraically as 
Pa{5p ~ Q,h — 0, t) ^ f". Simulating a lattice of linear 
size L ~ 4000 we obtain a = 0.41(1). From the scaling 
ansatz pa{L,t) — i~"^Ep^^(tL~^, 1), with z = v^^jv^. 
The best data collapse is obtained for z = 1.5(1). For 
more details see [T2| . 

DLG universality class ~ We are now able to address 
the question concerning which universality class the DLG 



belongs to. In table H] we compare the measured values 
of the DLG critical exponents with those of the Manna 
universality class [l3, showing that they are compatible 
with each other - for further details see [l3]- The shape 
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TABLE I: The measured critical exponents for the DLG and 
the corresponding critical exponents for the Manna universal- 
ity class in d = 2 [l^ . 



of the universal scaling functions is remarkably similar to 
those that can be found in the literature for the Manna 
universality class. This is seen e.g. by comparing our 
Fig. El with Fig. 4in [l3l, see also UM. From the scaling 
exponents and the scaling functions we find compelling 
evidence that the DLG belongs to the Manna universality 
class. 

Power spectrum at criticality ~ Finally we relate the 
exponent of the power spectrum at low density to the 
scaling properties of the DLG at the critical point pc 
of the APT[i3]. The scaling behavior of the correlation 
function Ca{r,t) = {pair,t)pa{0,0)) - (pa(0,0)2) at the 
critical density p = Pc determines the power law exponent 
in the scaling of the spectrum Saif) of the total number 
of active particles Na(t) in the lattice. Assume 



Cair,t) = X-^CaiXp^X-^'n^X-^^r) 



(5) 



In stationary directed percolation processes above criti- 
cality it has been observed [l^ that the correlation func- 
tion Ca(r, t) at t = first decays in space algebraically as 
j,-P/i'±^ until it saturates at a constant value at r > 
In the saturated regime the two sites become uncorre- 
lated so that this value is just equal to the squared sta- 
tionary density of active sites p^. One can use this fact 
to extract the scaling exponent for Sa{f), making use of 
the Wicner-Khinchin theorem: 

Ca{t) = {Na{t)Nam - {Na{Of) 

' I dPv [ dVp,(r,t)p„(r',0)\ + ... 
\Jv Jv I 



r 



so that: 

Saif) I Ca{t) e-^'^f' ^ /"'~^"" 



with p = 1 + i ^2 — j . With our estimates of the 
critical exponents for the DLG we find p = 1.82(6), 
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FIG. 5: DLG: Scaling behavior of the spectrum S{f) oc 
of the total number of particles in the box N{t) and of the 
spectrum Sa{f) of the total number of active particles near 
the critical density p = 0.24506 ~ Pc- The spectra have been 
multiplied by arbitrary factors to visualize the scaling expo- 
nent properly. Lattice size L — 1000. 



while with the best estimates of the Manna universal- 
ity class exponents [l^ one finds fi = 1.78(2). At low 
densities we expect that fluctuations in the total num- 
ber of particles N{t) in the lattice are triggered by active 
particles, therefore we expect the two spectra S{f) = 
^Jdt Nit) e-*2./t and SM) = ^ I dt Na{t) e'^^./* 
to show the same scaling behavior at p ~ pc, as it is 
confirmed by simulations (see Fig. [S]). Computing the 
power spectrum of the total number of particles in the 



box at p = 0.24506 ~ pc we observe a power spectrum 
compatible with the predictions, confirming that fluctua- 
tions in the total number of particles in the sub- volume V 
are determined by the scaling behavior of the correlation 
function at criticality. 

Summary - We have investigated the scaling prop- 
erties of the power spectrum S{f) of the total number 
of particles in a sub- volume of the lattice for the DLG: 
our simulations show that the DLG is not characterized 
by 1// fluctuations as had been previously observed, but 
present a much more variegated picture. At high densi- 
ties the power spectrum is the same as observed in a gas 
of random walkers, while at low densities the spectrum 
scales as S{f) ~ l//^'* and we have shown that this ex- 
ponent is determined by the decay of the density-density 
correlation function near criticality. 

We have shown that the deterministic lattice gas is in 
the same universality class as the Manna model[ia]. To 
our knowledge this is the first example of a completely 
deterministic and non-chaotic system in this universal- 
ity class. Perhaps this is not too surprising given that 
we found that a stochastic bulk noise is generated in 
the effective Langevin description as the system size is 
increased. From this point of view the DLG becomes 
similar to the other members of the Manna class, which 
all involved a stochastic element in their microscopic 
dynamics. Examples of models in the Manna class in- 
cludes the Oslo model[l^ and the Manna model[l^. In 
the Oslo model the local threshold for relaxation is up- 
dated stochastically at every relaxation and in the Manna 
model particles move to stochastically chosen neighbor 
sites. 
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